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Abstract. We carry out a systematic investigation on the low-energy electric dipole strength, 
which is often called pygmy dipole resonances (PDR), using the canonical-basis time-dependent 
Hartree-Fock-Bogoliubov (Cb-TDHFB) method. The Cb-TDHFB is a new method which is 
derived from TDHFB with an approximation analogous to the BCS theory that the pair 
potential is assumed to be diagonal in the time-dependent canonical basis. We apply the method 
to linear-response calculation for even-even nuclei. We report the neutron-number dependence 
of PDR in light (A < 70) and heavy isotopes (A > 100) around N = 82. 



1. Introduction 

Recent progress in radioactive isotope facilities and experimental techniques allows us to 
investigate properties of new exotic isotopes. In the neutron-rich nuclei, the neutron-skin and 
halo structures are known to be developed [I] . According to those exotic structures in the ground 
state, new types of elementary modes of excitation are expected. The low-lying electric dipole 
(El) modes, which are often called the pygmy dipole resonances (PDR), have been extensively 
studied as one of such characteristic modes of excitation in exotic nuclei [2]. However, the 
structure and the emergence mechanism of the PDR are still elusive at present. In this article, 
we perform a systematic investigation on the PDR to elucidate its characteristic features. 

In order to analyze and predict excited states in unstable nuclei, we use a time-dependent 
mean-field approach. The time-dependent Hartree-Fock (TDHF) theory has been extensively 
utilized to study nuclear many-body dynamics [3j. Recently it has been revisited with modern 
energy density functionals and more accurate description of nuclear properties has been achieved 
[H El El [U |9] . However, it neglects the pairing interactions in particle-particle and hole-hole 
channels, which are important for describing properties of open-shell heavy nuclei. It is well- 
known that the time-dependent Hartree-Fock-Bogoliubov (TDHFB) theory [10] properly takes 
into account the pairing correlation. Unfortunately, the practical calculation with the TDHFB 
are very limited [lit H2j . because it requires us to deal with the time evolution of "infinite" 
number of quasi-particle orbitals, in principle. In practice, the number of the quasi-particle 
orbitals is approximately equal to the dimension of the adopted one-body model space. 

Recently, we have developed the canonical-basis TDHFB (Cb-TDHFB) theory using a simple 
approximation to avoid to treat the "infinite" orbitals [13] . The Cb-TDHFB is a time-dependent 
version of the BCS approximation [J3] for the Hartree-Fock-Bogoliubov theory. Namely, we 



neglect off-diagonal elements of the pair potential in the time-dependent canonical-basis. In this 
paper, we use the Cb-TDHFB method to perform the linear-response calculation with the full 
Skyrme functional of the SkM* parameter set. We report the neutron number dependence of 
the low-lying El strength for light and heavy neutron-rich isotopes. 

The paper is organized as follows. In Sec.EJ we present the basic equations of the Cb-TDHFB 
method, their conservation lows and the numerical procedure of the present work. In Sec. [3j 
we present numerical results of the real-time calculations of the El mode for many isotopes. 
Finally, the summary is given in Sec. 21 



2. Cb-TDHFB equations 

2.1. Basic equations 

The Cb-TDHFB can be derived from the TDHFB equation on the assumption that the pair 
potential can be approximated to be diagonal in the canonical-basis. The readers should be 
referred to Ref. |13| for more details of the derivation. 

The TDHFB equation can be written in terms of the generalized density matrix R(t) and 
the HFB Hamiltonian U(t) as follows [TO] : 

i^R(t) = [H(t),R(t)}. (1) 

This is equivalent to the following equations for one-body density matrix p(t) and the pairing- 
tensor matrix n(t). 

i^p(t) = [h(t),p(t)] + K(t)A*(t)-A(t)K*(t), (2a) 

£ K (t) = h(t) K (t) + K(t)h*(t) + A(t){l-p*(t))-p(t)A(t). (2b) 
at 

Here, hit) and Ait) are the single-particle Hamiltonian and the pair potential, respectively. 
Introducing the time-dependent canonical states \(f>k(t)) and \(j>k(t)}, we express the TDHFB 
state in the canonical (BCS) form as 

i*(*)> = n + «*(*)4w4(*)} i°>- ^ 

Here, it should be noted that the pair of states, k and k, are not necessarily relate to each other 
by the time-reversal, \4>k(t)) 7^ T\(j)k(t)). The diagonal approximation of the pair potential lead 
Eqs. (|2aj) and (f2b|) to the following equations: 

iftlMt)) = (Ht) - Vk(t))\<t>k(t)), iftlfait)) = (Ht) ~ Tk(t))\<Pk(t)), (4a) 
ij t Pk(t) = Kk(t)A* k (t) - «Ut)A fc (t), (4b) 
*^Kfc(t) = (Tfc(t) + Visit)) K k {t) + A fc (t) (2p k (t) - 1) . (4c) 

These basic equations determine the time evolution of the canonical states, and \(j>k(t)), 

their occupation pk(t) = \vk(t)\ 2 , and pair probabilities Kk(t) = Uk(t)vk(t). The time-dependent 
gap Afc(t) is calculated as 

M*) = -£«j(*)«fciMr, ( 5 ) 

l>0 



using the anti-symmetrized two-body interaction v |13] . The real function of time nt (t) is defined 

as 

Vk(t) = (Mt)\h(t)\Mt)) + i C^r\<t>k(t)\ . (6) 

This is an arbitrary gauge parameter and depends on the phase of the canonical states. Note 
that the symbols, (p, n, A), indicate matrices with two indices in Eqs. (I2ap and (|2bj) . but vectors 
with a single index in Eqs. (|4ap . (|4bj) . and (@c|). Similar equations can be found in Ref. [15] for 
a very schematic pairing energy functional. 



2.2. Pairing energy functional and cut-off function 

We adopt a Skyrme functional with the SkM* parameter set for the particle- hole channels. For 
the pairing energy functional, we adopt a simple functional of a form 

£pair(t) =~J2 Gkl4(t)Kl{t) =-J2 Kfc(t)Afc(t), A fc (t) = GkMt), (7) 
k,l>0 k>0 Z>0 

with a following special gauge condition [13], 

= £fc(t) = (Mt)\h(t)\Mt)), Vk(t) = = <<^(*)IM*)I<^(*)>- (8) 

Here, G^/ = Gf{e^)f{e®) and the constant G is determined by the smoothed pairing method 
[13] . The cut-off function /(e?) depends on the single-particle energy of the canonical state k 
at the HF+BCS ground state. The cut-off function /(e) is written as 

9(e c -s), (9) 
with the cut-off energies 

e c = A + 5.0 MeV, e c = e c + 2.3 MeV, (10) 

where A is the average of the highest occupied level and the lowest unoccupied level in the HF 
state. Here, the cut-off parameter e c is necessary to prevent occupation of spatially unlocalized 
single-particle states, known as the problem of unphysical gas near the drip line. For neutrons, 
if e c becomes positive, we replace it by zero [T6] . 



1 + exp 



0.5 MeV 



2. 3. Properties of the Cb- TDHFB equations 

The Cb-TDHFB equations possess the following properties (13] : 

(i) Conservation law 

(a) Conservation of orthonormal property of the canonical states 

(b) Conservation of average particle number 

(c) Conservation of average total energy 

(ii) The stationary solution corresponds to the HF+BCS solution, 
(hi) In the limit of A = 0, they are equivalent to TDHF. 

(iv) In the small-amplitude limit, if the ground state is in the normal phase, the equations 
are identical to the particle-hole, particle-particle, and hole-hole RPA with the BCS 
approximation. 



2-4- Numerical details 

For numerical calculations, we extended the computational program of the TDHF in the 
three-dimensional (3D) coordinate-space representation [2] to include the BCS-type pairing 
correlations. The ground state is first constructed by the HF+BCS calculation. Then, we 
add a weak impulse isovector dipole field V E i(t) to the ground state. V E i(t) is written as 

. . „, . I (Ne/A)n for protons . . 

V E i(t) = - V F E1 8(t), F E i = C ' .! P , (11) 

\—{Ze/A)ri tor neutrons 

where i = (x, y, z). We solve the Cb-TDHFB equations in real time and real space, then calculate 
the linear-response of the nucleus. We calculate the time evolution of the expectation value of 
Fex under the external field V E i(t) with rj <^ 1 and obtain the strength function S(E1;E) by 
Fourier transformation using an exponential smoothing with T = 1 MeV: Vei(^) — t VEi(t)e~ Tt l 2 . 

In this work, we use the 3D Cartesian coordinate for the canonical states, (j)i(r,o~;t) = 
(r , cr\(f>i(t)) with a = ±1/2. For light isotopes (A < 70), the coordinate space is discretized in 
the square mesh of Ax = Ay = Az = 0.8 fm in a sphere with radius of 12 fm. For heavy 
isotopes (A > 100), we choose the square mesh of Ax = Ay = Az = 1.0 fm in a sphere with 
radius of 15 fm. 



3. Neutron-number dependence of the low-lying El strength 

In this section, we discuss a systematic trend of the low-lying El strength as a function of 
neutron number. To quantify the strength of PDR, the following ratio of the low-lying El 
strength is used: 



mi(E c 



S(El;E)dE 



x 100 



mi " / E X S(El;E)dE 
where E c is a cut-off energy. We take E c =10 MeV for all calculations. 



(12) 
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Figure 1. Neutron number dependence of 
m\{E c )/mi defined in Eq.((T2J| for O, Ne and Mg 
isotopes 
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Figure 2. Same as Fig. 1 but for S, Ar and Ca 

isotopes. 



3.1. Light to medium nuclei (A < 70) 

For light nuclei with Z < 40, the systematic analysis has been carried out in Ref. |17j . 
The calculation was based on the finite amplitude method \18\ [T9] with the Skyrme energy 
functional of SkM*, however, the pairing correlation was neglected. Thus, it would be desirable 
to investigate the effect of the pairing correlations on the low-energy El strength. 
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Figure 3. Same as Fig. 1 but for Sr, Zr, Mo, Ru, Pd, Cd, Sn, Te and Xe isotopes (N = 78 - 
86) are shown as a function of neutron number. All results are obtained by using Cb-TDHFB. 

Figure 1 shows the results for O, Ne and Mg isotopes as a function of neutron number and 
Fig. 2 shows those for S, Ar and Ca isotopes. Solid lines with filled symbols indicate the results 
of Cb-TDHFB calculation and dashed lines with open symbols are those of Hartree-Fock plus 
random phase approximation (HF+RPA) without the pairing correlation [T7J. In both cases, 
we see that the nuclei with A" =8 - 14 have the fractions of low-lying El strength less than 
1.0%, and a sudden jump of the fraction at N = 14 — >• 16 on each isotopic chain in Fig. 1. 
The neutron number A^ = 16 corresponds to the occupation of the S1/2 orbit. Figure 1 shows a 
neutron shell effect on low-lying El strength. This suggests that the position of the Fermi level 
plays an important role for the emergence of low-lying El strength. 

A similar behavior is also seen in Fig. 2, however, is slightly changed by the pairing 
correlations. Namely, the nuclei with A" =16 - 28 have small fractions in HF+RPA, and the 
sudden jump from A^ = 28 to 30 is clearly observed. In Ca isotopes, the fraction has almost 
constant values at A^ > 36. The neutron number A^ = 30 corresponds to the occupation of the 
P3/2 orbit, and A^ = 36 to the orbit. The occupation of the orbits with small orbital angular 
momentum is important for the emergence of low-lying El strength. 

The effects of the pairing correlation is mostly seen around N = 28 in Fig. 2. The low- 
lying El strength increases smoothly around A" = 28 in S and Ar isotopes in the Cb-TDHFB 
calculation, while we see sudden increases from A^ = 28 to 30 in the calculation without the 
pairing (HF+RPA). This smooth evolution of the low- lying strength is caused by the fractional 
occupation probabilities of the single-particle orbits, due to the pairing correlation. For lighter 
nuclei in Fig. 1, the pairing effect is seen in Ne and Mg isotopes. This is mostly due to different 
shapes in the ground state. Actually, for AT = 16 — 20, the calculated ground states in the 
HF+BCS are spherical, but those in the HF are deformed. 

3.2. Heavy nuclei {A > 100) 

Figure 3 shows the results for Sr, Zr, Mo, Ru, Pd, Cd, Sn, Te and Xe isotopes with N = 78 - 86, 
calculated with Cb-TDHFB. Again, we can clearly see a sudden increase of the PDR ratio from 
A^ =82 to 84, for all isotopic chains. The kink position corresponds to the well-known neutron 
magic number N = 82. Even in the case of heavy isotopes, the neutron shell effects play a very 



important role in the low-lying El strength. However, there are constant but significant PDR 
fractions for N =78 - 82. Although 132 ~ 136 Xe are stable nuclei, they have 1.5% of low-lying El 
strength. This fraction even increases as the proton number decreases. We suppose that these 
constant fractions have a different origin from the one observed in the light neutron-rich nuclei. 
This is currently under investigation. 



4. Summary 

We presented an approximate and feasible approach to the TDHFB; the canonical-basis TDHFB 
method. We calculated the El strength distribution in neutron-rich isotopes of light and heavy 
nuclei, using the real-time real-space method. The low-lying El strength has a characteristic 
neutron number dependence. In light neutron-rich isotopes with Z = 8 - 20, the fraction of 
the energy- weighted El strength below 10 MeV increases at N = 16 in O, Ne and Mg isotopes 
and at N = 30 in S, Ar and Ca isotopes. These numbers indicate that the occupation of the 
single-particle states with small orbital angular momentum is very important for the emergence 
of the low-lying El strength. The heavy nuclei with neutron excess around A = 82 have been 
also investigated. There is a similar jump in the PDR fraction at N = 82 — > 84. However, a 
significant amount of the PDR fraction is also observed in nuclei before the jump with N < 82. 
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